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ABSTRACT 


This  paper  is  primarily  concerned  with  selecting  a  subset 
of  k  populations  such  that  the  probability  is  at  least 
P  that  the  selected  subset  includes  the  population  with 
the  largest  (smallest)  quantile  of  a  given  order  a 
^  a  1)  •  particular  a  procedure  is  proposed  and 
studied  which  is  valid  for  any  family  of  distributions 
with  increasing  failure  rate  on  the  average  (IFRA) .  It 
is  compared,  asymptotically,  with  a  distribution-free 
procedure  proposed  by  Rizvi  and  Sobel. 


SELECTION  PROCEDURES  FOR  RESTRICTED 
FAMILIES  OF  PROBABILITY  DISTRIBUTIONS+ 


by 

Richard  E.  Barlow  and  Shanti  S.  Gupta 

1.  INTRODUCTION  AND  SUMMARY 

Let  ni,n2,  . . . ,  be  k  populations.  The  random  variable 
associated  with  11^  has  a  continuous  distribution  F^,i  ■  1,2,  k  .  We 

are  primarily  interested  in  selecting  a  subset  such  that  the  probability  is  at 
least  P*  that  the  selected  subset  includes  the  population  with  the  largest 
(smallest)  quantile  of  a  given  order  a(0  <  a  <  1)  .  We  assume  each  F^  has 
a  unique  a-quantile,  ^  .  Let  F^(x)  •=  F^j  denote  the  cumulative  distri¬ 

bution  function  of  the  population  with  the  ith  smallest  a-quantile.  In  the 
following,  we  consider  families  of  distributions  ordered  in  a  certain  sense  with 
respect  to  a  specified  continuous  distribution  G  and  propose  and  study  a 
selection  procedure  which  is  different  from  the  nonparametric  procedure  of  Rizvi 
and  Sobel  (1967).  We  assume 


(a)  F^Cx)  l  Fjkj(x)  i  *  1,2,  ...,  k  and  all  x  , 

(b)  3  a  continuous  distribution  G  3  F^j  *<  G,  Vi  *  1,2,  ...,  k  , 


where  •(  denotes  a  partial  ordering  relation  on  the  space  of  distributions. 


This  research  was  partly  supported  by  the  Office  of  Naval  Research  contract 
N0NR-1100(26)  at  Purdue  University  and  Office  of  Naval  Research  contract 
NONR-3656(18)  at  the  University  of  California  at  Berkeley.  Reproduction  in 
whole  or  in  part  is  permitted  for  any  purposes  of  the  United  States  Government. 
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A  relation  on  the  space  of  distribution  is  a  partial  ordering  if 

F  ^  F  V  distributions  F 
F  «  G,  G  <  H  implies  F  <  H  . 

Note  that  F  G  and  G  F  do  not  necessarily  imply  F  E  G  . 

Various  special  cases  in  audition  to  stochastic  ordering  are: 

F  ■<  G  iff  F (0)  -  G (0)  -  0  and  is  (i) 

*  x 

nondecreasing  in  x  >  0  on  the  support  of  F  . 

F  -K  G  iff  G  ^F(x)  is  convex  on  the  support  of  F  .  (ii) 

c 

F  <  G  iff  F (0)  -  G(0)  -  Jj  and  is  (iii) 

r 

increasing  (decreasing)  for  x  positive  (negative)  on  the  support  of  F  . 

* 

F  G  iff  F(0)  -  G(Q)  ■  *s  and  G  *F  is  concave-convex  (iv) 

s 

about  the  origin,  on  the  support  of  F  . 

If  G(x)  ■  1  -  e  X  for  x  >  0  ,  then  (i)  defines  the  class  of  IFRA 
distributions  studied  by  Birnbaum,  Esary  and  Marshall  (1966)  while  (ii)  defines 
the  class  of  IFR  distributions  studied  by  Barlow,  Marshall  and  Proschan  (1963). 

For  any  distribution  G,  F  4  G  iff  F(x)  crosses  G(0x)  at  most  once  and  from 

* 

below  if  at  all  as  a  function  of  x  for  all  9  >  0  .  If  G(x)  «  1  -  exp(-x  ) 

for  x  >  0  and  A  >  0  ,  then  F  -K  G  implies  that  F  is  "sharper"  than  the 

* 

family  of  Welbull  distributions  with  shape  parameter  A  .  Implications  of 
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orderings  defined  by  (iii)  were  studied  by  Lawrence  (1966).  Van  Zwet  (1964) 

studied  orderings  defined  by  both  (ii)  and  (iv).  Clearly  ordering  implies 

c 

ordering  and  ordering  implies  K  ordering. 

*  s  r 

If  =  (X^jX^*  ***»  ^in^  *s  t*'e  °bserved  sample  from  the  ith  population, 
then  we  restrict  ourselves  to  the  class  of  statistics  *=  T(X^)  thac  preserve 
both  ordering  relations  (a)  and  (b) ,  i.e., 

P„  (T(X)  <  x)  >  P  (T(X)  <  x}  (a') 

Fu)  "  ‘  "  m 

for  all  x  and  i  =  1,2,  ...»  k  . 


FY(X1)  ^gt(y) 


(b') 


i  *  1,2,  ....  k  where  Y  *  (Y.,Y_,  ....  Y  )  is  a  random  sample  from  G  . 

-it  n 

In  Section  2  of  this  paper,  we  propose  and  study  procedures  R  (R‘)  for 
selecting  the  population  with  the  largest  (smallest)  a-quantile  for  distributions 
which  are  •*(  ordered  with  respect  to  a  specified  distribution  G  .  The  infimum 


* 

of  the  probability  of  a  correct  selection  is  obtained  In  Theorem  2.1  and  asymptotic 

evaluation  is  given  in  Theorem  2.2.  Section  3  deals  with  quantile  selection 

procedures  for  the  class  of  IFRA  distributions.  In  Section  4,  we  study  the 

efficiency  of  procedure  R  with  respect  to  a  procedure  studied  by  Rizvi  and  Sobel 

(1967)  under  scale  type  slippage  configurations.  Asymptotic  relative  efficiency 

of  R  with  respect  to  a  selection  procedure  for  the  gamma  populations  proposed 

by  Gupta  (1963)  is  also  investigated.  Section  5  deals  with  selection  procedures 

for  the  median  for  distributions  that  are  ordered  with  respect  to  a  specified 

r 

G  .  In  Section  6  we  propose  a  selection  procedure  with  respect  to  the  sample 
means  for  distributions  that  are  ^  ordered  with  respect  to  C(x)  c  1  -  e  *  . 


c 

Appl  1  cation  to  the  selection  of  gamma  populations  is  also  given  In  Section  6. 
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2_. _ QU ANTT1.K  SELECTION  RIHLFS  FOR  DISTRI HUTIONS  -<  ORDKRKD  WITH  RKSPKCT  TO  G 

* 

We  are  given  a  sample  of  size  n  from  each  of  the  k  populations 

ni,i  =  1,2,  ...,  k  .  The  distributions  and  the  specified  distribution 

G  satisfy  the  assumptions  (a)  and  (b)  of  Section  1.  The  distributions 

are,  otherwise,  unspecified.  Of  course,  the  correct  pairing  of  the  unordered 

and  ordered  F^'s  Is  not  known.  We  denote  the  k-tuples  *k^  ^y 

0  .  Let  T  denote  the  jth  order  statistic  from  F  where  j  <  na  <  j  +  1  . 

J » i  1 

Clearly,  T  -*■  (,  .  ,  the  a-quantile  as  n  ■+  00  and  o  .  The  rule  we 

1 ,  i  a ,  i  M  n 

J  a.s. 

propose  for  selecting  the  population  with  the  largest  a-quantile  is: 

R  :  Select  population  iff 

T  >  c  max  T  ,  j  <  no  <  n  +  1  ,  (2.1) 

J*1  ~  1  <  r  <  k  ^’r 

where  c  *  c(k,  P*,  n,  j)  is  some  number  between  0  and  1  which  is  determined 
so  as  to  satisfy  the  probability  requirement 

inf  P{CS | R>  -  P*  ,  (2.2) 

n 


where  CS  .stands  for  a  correct  selection,  i.e.,  the  selection  of  any  subset 
which  contains  the  population  with  distribution  F^j  •  Before  discussing 

the  main  theorem  concerned  with  the  evaluation  of  P { C S | F. }  ,  we  present  a  known 


result  for  order  statistics.  Let  H  (x)  be  the  cdf  of  the  jth  order  statistic 

J » 1 

from  F^j  and  let  (x)  be  the  cdf  of  the  jth  order  statistic  from  G  .  Let 
us  define 


_ nj _ 

( j-1) ! (n-j )  ! 


f  ul-1  (1  -  u)n-J 


du 
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so  that 


'  BJ,n<F[i]°t)>  5  BJ,nF[ilW  •  (2'3) 

Since 

'  Cj1  H  ^(x)  -  t BjG ] _1  B  F^j  (x)  -  G"1  Fu)(x)  ,  (2.4) 

we  see  that  order  statistics  preserve  each  of  the  partial  ordering  relations 
(i)  -  (iv) .  For  additional  applications  of  (2.4)  sec  van  Zwet  (1964). 

We  now  state  and  prove  a  theorem  which  enables  us  to  compute  the  constant 
c  which  defines  the  procedure  R  . 

Theorem  2.1: 

If  F(i](0)  -  G (0)  -  0,  Fu](x)  >  F(k](x),  x  >  0,  i  -  1,2 . k  and 

F[ki  ■  then 

Inf  P(CS|R)  -  j  [G  (a))k~1  dC,(x)  .  (2.5) 

a  o  j  c  J 


Proof : 

Note  that 


p{cs|r} 


dHj.k(x) 


o 


dHJ  kW  . 
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00 

-  /  [Gj  (^f)  ]k_1  dGj  (x)  , 

provided  c  is  between  0  and  1.  | | 

Remark  1: 

The  constant  c  ■  c(k,  P*,  n,  j)  which  defines  the  selection  procedure  R 


is  determined  by 
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oo 

/  lGJ(*)]k'1  dG  j  (x)  =  P*.  (£  <  P*  <  1) 


(2.9) 


These  constants  arc  tabulated  for  G(x)  =  1  -  e  for  selected  values  of  n,  k,  j 
and  P*  in  the  first  set  of  tables  in  the  companion  paper  by  Barlow,  Gupta 
and  Panchapakesan  (1967).  Clearly,  c  is  independent  of  scale. 


Remark  2: 


-(-)X 

V 


If  Gj (x)  *’ 1  -  e  ,  for  x  >  0  and  9,  X  >  0  ,  then  for  j  *  1  ,  the 
values  of  c  are  independent  of  n  .  This  can  be  seen  from  the  fact  that  the 
distribution  of  the  smallest  order  statistic  involves  n  only  as  a  scale  parameter 
and  that  the  selection  procedure  (2.1)  is  scale  invariant. 


Remark  3: 


It  should  be  pointed  out  that  Theorem  2.1  requires  only  F^j  G  ; 
however,  to  apply  the  procedure  R  ,  we  assume  that  F[i]  ^  G*  Vi  * 

Now  we  discuss  the  asymptotic  evaluation  of  the  probability  of  a  correct 
selection.  We  state  and  prove  the  following  theorem. 


Theorem  2.2: 


If  F[it]  (*)  G,  F^j  (G)  has  a  differentiable  density  f^j  (?)  1°  a 

neighborhood  of  the  a-quantile  £  (n  )  and  fr.  .(£  )  i  0  (g(n  )  +  0)  ,  then  in 

o  a  l  k  J  a  a 

our  previous  notation  (see  Theorem  2.1) 


.  (  *1.k  )  -f  k-l/x  ,  W*  \ 

„  max  X,  ^  C  “  J  *  \7  +  - 7= - J 

J’k(l<r<k  j'r  )  ' 


d$  (x) 


d$  (x) 


(2.10) 


V 
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as  n  -*•  °°  ,  where  a  =  1  -  a  and  4>(*)  is  the  c.d.f.  of  the  standard  normal 
variate. 

Proof : 


P)X,  ,  >  c  max  X 


Ik- 

J  ’  1< 


r<k-l 


j.r 


c-1 


max  (X,  -C  )  +  —  5 

,  e  3,r  *  c  * 
/d  8  /a  a/(/n  fri1(£  )) 


lk]'V 


(2.11) 


,  (1~c)  hki  <« 


since 


•  l 


d<t(x)  , 


To  prove  the  second  part  of  (2.10),  note  that  F,, ,  G  implies 

l  ^  J  uu 

-1 

G  F^j(x)  -  x  changes  sign  at  most  once  and  from  -  to  +  ,  if  at  all.  Since 

Fr.  ,  G  is  invariant  under  scale  changes,  we  can  assume  n  ■  £  so  that 
[k]  A  6  ’  a  a 

F(£^)  £  G(^a)  •  Either  F  »G  in  a  neighborhood  of  or  the  slope  of  the 

tangent  line  to  F  at  is  greater  than  the  slope  of  the  tangent  line  to 

G  at  .  In  either  case  f^j  (S^)  >  g(£Q)  and  in  general  5^  f^j  (S^)  > 

\  g(na)  * 


Remark  4: 


Setting 


\c  c/r^  / 


)  -  p* 


(2.12) 
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The  values  of  d  are  tabulated  for  selected  values  of  k,  P*,  n  and  j  in  the 
companion  paper  by  Barlow,  Gupta  and  Panchapnkesan  (1967). 

The  rules  R  and  R'  select  nonempty  subsets.  The  size  of  the  selected 

subset  is  a  random  variable  which  takes  values  1,2 . k  .  The  expected  size 

of  the  selected  subset  is  a  common  measure  of  the  efficiency  of  the  procedure 
(Gupta  (1963b)).  However,  it  is  difficult  in  our  more  general  framework  to  set 
meaningful  bounds  on  the  expected  size  without  further  assumptions.  If  we 

assume,  in  addition,  that  there  exists  G+  such  that  G+"Kf.  .  for  all  i  , 

* 

then  we  can  obtain  an  upper  bound  on  the  probability  of  including  the  "worst" 
population  in  the  selected  subset.  We  consider  this  in  more  detail  later  for 
IFRA  distributions. 

If  we  assume  that  F^^x)  stochastically  increasing  with  respect  to 
i  ,  then 

Ptselect  |R)  >  P{select  n^j|R}  if  i>j  .  (2.17) 

The  proof  is  similar  to  the  one  given  in  Gupta  (1967).  A  result  similar  to 
(2.17)  is  true  for  R^  . 
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3_. QUANTILE  SELECTION  PROCEDURES  FOR  THE  CLASS  OF  I  FRA  DISTRIBUTIONS 

If  F  *<(  G  where  G(x)  -  1  -  e  X  for  x  >  0  ,  then  F  Is  an  IFRA  distri- 
* 

butlon.  The  problem  of  selecting  the  best  one  of  several  IFRA  populations  has 
been  considered  by  J.  K.  Patel  (1967).  He  was  interested  in  selecting  that 
population  with  the  smallest  failure  rate  at  a  prescribed  time  T  .  His  decision 
rule  depends  only  on  the  number  of  observed  failures  in  [0,  T]  for  each 
population  and  not  on  the  times  at  which  failure  occurred.  In  this  respect  our 
procedure  will  utilize  more  information,  though  of  course  we  are  selecting  with 
respect  to  quantiles  rather  than  failure  rate. 

We  show  how  to  obtain  c  values  for  the  Weibull  distribution  with  shape 

A 

parameter  \  >  0  .  We  remark  that  the  class  of  distributions  F  ,  such  that 

X 

x 

F  ^  G  where  G.  (x)  *  1  -  e  ®  for  x  >  0  and  6,  X  >  0  is  the  smallest 

i  A  A  " 
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or 


If,  in  addition  to  the  assumptions  of  Section  2  (see  Theorem  2.1),  we 

assume  that  a)  F^^x)  2  F^j(x)  ^or  x  t  *  =  1*^ . ^  an<* » 

b)  G  ^  F,  ,  for  all  i  =  1, 2,  , . . ,  k,  A  >  1  then  we  can  obtain  an  upper  bound 
on  the  probability  of  selecting  the  "worst"  population,  i.e.. 


.  P{ Selecting  n^^R}  <  /  {Gj(-^)lk“1  dGj  (x) 


(3.2) 


where  c  is  chosen  so  that 

00 

P{CS|r}  >  /  [Gj  (^)  l^"1  dGj  (x)  -  P*  . 

Clearly,  the  upper  bound  is  an  increasing  function  of  X  for  X  >  1  . 

Now  we  describe  the  relative  performance  for  small  sample  size  (n  *  15) 
for  procedures  R  and  R^  (the  Rizvi-Sobel  procedure  -  see  (4.8))  using  Monte 
Carlo  techniques.  For  this  purpose  we  chose  gamma  and  Weibull  distributions 
with  densities 


gamma. 


i  -  1,2 


Weibull, 


i  ■  1,2 


Based  on  5,000  simulations,  we  computed  P ( C S J R } ,  E (S | R)  and  P{CS|R^}  , 

E(S t Rx )  .  (  E(S|R)  is  the  expected  size  of  the  selected  subset  using  procedure  R 

These  values  are  given  below: 


I  v 
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Monte  Carlo  Conparisons^of  R  and 


P*  =  .90,  k  =  2,  n  =  15 


Gamma 

Gamma 

Weibull 

r  -  1 

r  =  5 

r  =  2 

\  =  1.  S2  "  ^ 

0-L  =  2,  92  =  3 

9!  =  1.  92  "  2 

r{cs I k} 

.993 

1.000 

.985 

e(s!r) 

1.47 

1.88 

1.85 

p{cs ! R  } 

.997 

1.000 

•  939 

f.(s,r1  ) 

1.64 

1.96 

1.76 

Computations  performed  by  David  Stanford  of  the  Operations  Research  Center, 
University  of  California,  Berkeley. 
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4.  EFFICIENCY  OF  PROCEDURE  R  UNDER  SLIPPAGE  CONFIGURATIONS 

We  consider  slippage  configurations  F^j(x)  =  F(— ),  i  =  1»2,  ....  k  -  1 
and  F[^(xJ  =  F(x),  0  <  6  <  1  .  We  obtain  asymptotic  expressions  for  the 
probability  of  a  correct  selection  and  the  expected  size  of  the  selected  subset 
for  procedure  R  and  for  two  other  procedures. 

Using  our  previous  notation 


P{CS|r}  ■=  P(T  ,  >  c  max  T  ) 
I  "  l<r<k  j,rj 


j.k 


X,  ,  >  c  max  5  X. 
Jlk  -  l<r<k 


(4.1) 


-  PG  Yi  k  -  Cfi  “aX  Yi  r 
Gj  I  J’k  l<r<k  J’r) 


where  Y  ,  r  *  1,2,  k  are  i.i.d.  with  the  cdf  G  (y)  .  From  (2.11),  we 

J  • r  J 

obtain 


"  L.  i  /  v  (l-c5K  f(£  )/n\ 

P{ CS  | R }  =  /  *k_1  +  - - - - - I  dt(x) 

\  c6/a  5  / 


(4.2) 


Note  that  the  probability  of  a  correct  selection  is  a  monotone  decreasing 
function  of  6  .  For  the  slippage  configuration 


E (S | R)  -  P(CS | R)  +  (k  -  1)  PIT  .  >  c  max  T  .  j 

(  '  i*i 


(4.3) 


PlT  .  >  c  max  T 
"  1*1 


-  P 


It  -6£  >  c  max/  max  (T  -6E  ),T.  .< +*-«€„  \+c6£  -6E  J  (4.4) 

!  J’  \2<i<k-i  j>i  °  j*k  a  a  a)  a  aj 


-  L  t(r  K  (1-  -)/T,\  /  r  f (r  )(l  -  -)/Ti\ 

«  f  - £_J.>k-2f£._H - 1 - - — - - ) 

-<*>  \  c  /a  u  /  \  /ii  5  ' 


d : (x)  . 
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Sotting  k  =  2  ,  wo  have 


E(S I R)  -  P{CS|R} «  ♦  I- 


f«„)  ttt(i  -~)/i 

/m  -Ji  *  c-)2 

c 


(A. 5) 


Setting  the  right-hand  side  of  (4.5)  equal  to  e,  -  f  (CQ)  ^1  - 

/a  a-.  /]  +  ' \2  4>  (e)  ,  where  we  have  put  c  *  c  .  Now  using  c  « 

V  \Cn/ 


i-  ^ 


(from  (2.14)) 


-  f(C 


a)  ^1-6-^  +  62(i  +  +  2C' 


(4.6) 


from  which,  keeping  terms  of  order  </n  ,  we  obtain 


[  -  r1(e)(l  +  62)15  +  /2C6£  f(£  )]2 


(e)  * 


2  2  2 

r  ra  mi  -  6)^ 

a  a 


(4.7) 


Comparison  with  Rizvl  -  Sobel  Procedure 

Rizvi  and  Sobel  (1967)  propose  and  investigate  a  distribution-free  procedure, 
,  for  the  quantile  selection  problem. 

R^  :  Select  population  iff 


T.  .  >  max  T, 

-1’1  -  l<r<k  3'a’r 


(4.8) 


where  a  is  the  smallest  integer  with  1  <  a  <  j  -  1  for  which 


inf  P{CS ! R  }  >  P* 

n 


(4.9) 


is  satisfied. 


A  disadvantage  of  this  procedure  is  that  for  any  given  a  and  k  a  value 
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of  a  <  j  -  1  may  not  exist  for  some  pairs  (n,  P*)  .  However  if  P*  is 
chosen  not  greater  than  some  function  P^(n,  a,  k)  where  1/k  <  P^  <  1  , 
then  a  value  of  a  <  j  -  1  does  exist  that  satisfies  (4.9).  Rizvi  and  Sobel 
compare  the  efficiency  of  this  procedure  relative  to  several  competing  procedures 
under  translation  configurations. 

We  discuss  the  asymptotic  probability  of  a  correct  selection  using  their 
procedure  under  the  slippage  configuration. 


p(cs!r  1 


P  Tj  t  >  mox  T. 

I  3-k  •  l<r<k-l  3' 


a,  r 


a 


J  -6C  +(6-l)£  \ 

—  >  max  (-^ — 1 - - - J  6/n  f(6  ) 

l<r<k-l  \  6/a  a  /  a 


(4.10) 


-/ 


*  4. 

\«  7-» 

_oo  \  /a  a 


v  j 


s/rs 


di>  (x)  . 


The  derivation  above  uses  the  fact  that 


T.  -56 

-taa: — a  /-«  I  f(5 

/a  ii  0  law 


l) 


where  y//n  =  a/(n  +  1)  .  (See  Lemma  2  of  Rizvi  and  Sobel  (19671) 
Simiar lv 

E(S[R1)  -  P{ CS I Rx } ~ 


(1-6)6  /5  f (6  ) 

a  a 


d>(x) 


(4.11) 


Settinc  k  =  2  ,  we  have  from  (4.11) 
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E(S|R,)  -  p(cs|V  *  J— :{)  wjv*  \ 

'  /cTS  n/T+  62  / 


(4.12) 


Equating  the  right-hand  side  of  (4.12)  to  e  ,  we  obtain 


q\4+62  t’~1(c)-vl2 


n  (e)  %  -I--0  aVI+6 — * _ lfJblJ 

1  (l-i)2[^  f(^))2 


(4.13) 


For  the  slippage  configuration  above  we  define  the  asymptotic  relative 
efficiency  A  R  ECR.R^fi)  of  R  relative  to  Rj  to  be  the  limit  a.  c  *  0 


of  the  ratio  of  n  (e) 


t0  Ve> 


"r 

A  R  E(R,R  ;  6)  -  lim  -1-;  ■ 

1  eiO  Ve) 


VoTl  ^1 »6 2  -  y]2 _ 

$  1(e)\/(l +62)  +  sfi&O,  f  (£  )12 

a  a  J 


(4.14) 


Using  the  fact  that  F  -<  G  implies 


*a  f(5a)  -  na  *(na>  ’ 


we  see  that 


pf{cs|r1>  >  pg(cs [r1j 


(4.15) 


ef(s|r1)  -  pf(cs J Rl }  <  Ec ( s | r ^ )  -  p^fcr.jR^ 


(4.16) 
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where  both  (4.15)  and  (4.16)  are  asymptotically  true  as  n  -*•  •  for  the 
slippage  configuration. 


Comparison  with  Gupta  Procedure 


Gupta  (1963)  gave  a  selection  procedure  for  gamma  populations  with 


densities 


exp(-  f-)  if-)*"1,  x  >  0,  01  >  0,  i  -  1,2 . k 


This  p-ocedure,  ,  based  on  the  means  of  sample  size  n  from  each  of  the  k 
populations  is: 


R2:  Select  the  population  corresponding  to  the  observed  mean  x^  iff 


x.  >  b  max  x. 
1  ”  l<j<k  J 


(4.17) 


where  b  is  the  largest  constant  (0  <  b  <  1)  chosen  so  that  PfCS^}  >  P*  . 
Letting  v  -  2n  r  ,  it  is  shown  that  logg  b  ss  -d  where  d  is  independent 


of  n  and  satisfies 


J  ^  (x  +  d)  d$(x)  ■  P* 


(4.18) 


Assume  that  the  ranked  9^'a  have  the  slippage  configuration  9[i]  '  66(k] 


0  <  6  <  1,  i  •  1,2,  ...,  k-  1  .  Then 


E(S|R„)  -  PfCSlRj  a  (k  -  1)  f  *k'2/x-  — $/x  -  d?(x) 

22  -•  \  m^Dj  \  m^D )  .. 


(4.19) 


so  that  for  k  ■  2 
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e(s|r2)  -  p(cs|r2) 


{ 


ZJQ£  (b/6) 


) 


#0 


(4.20) 


Setting  the  right-hand  side  of  (4.20)  equal  to  e  and  solving  for 

n  =  n  (e) 

R2 


„  (c)  a12£hsl  -  /2]2 

2r(log  6)2  ' 


(4.21) 


\  (e) 

A  R  E(R,R  •  6)  %  lim  — ■  - 
2  CIO  V0 


2(l-6)2[Ca  f(Ca)}2 
r (log  6 ] 2  a  a(l+52) 


(4.22) 


for  r  >  1  , 


A  R  E(R,R2;6)  >  j-U~a)2  U~a)2  f-log  (l-oPT 

r(log  6)2  a  a (1+52) 


f(C  )] 

Also  lim  A  R  E(R,R  ;<5)  -  —2 - -2 - 

<5U  1  r  a  a 


(4.23) 


so 


that  lim  A  R  EIR.R^S)  >  11’”)1 \  (UttIng  r 


1) 


For  a  -  4  and  r  -  1,  lim  A  R  E(R, R  ; 6)  -  (log  2]2  -  .493  . 

6+  1 
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5 •  SF1  KCTION  WITH  RKSPLCT  TO  THE  MEDIAN  FOR  DISTRIBUTIONS  r-ORDF.KLD  WITH 
RKSPECT  TO  A  SPECIFIED  DISTRIBUTION  G 

We  consider  selection  procedures  with  respect  to  the  median  for  distributions 

which  have  lighter  tails  than  a  specified  distribution  G  .  We  say  that 

F.  has  lighter  tails  than  G  if  F.  centered  at  its  median,  A.  ,  is 
1  i  i 

4  -  ordered  w'ith  respect  to  G  (G(0)  =  h)  and  ~  F .  (x  +  A.)  I  >  ~  G(x)|  _  . 

r  *  dx  i  i  'x=0  -  dx  'x=0 

Here  we  are  following  an  ordering  proposed  by  Doksum  (1967).  (See  van  Zwet  (1964), 
Lawrence  (1966).) 

We  wish  to  select  a  subset  of  the  k  populations  containing  the  population 
with  the  largest  median  A^j  .  The  selection  rule,  we  propose,  is  in  terms  of 

the  sample  medians.  We  use  the  same  notation  as  in  Section  2.  The  rule,  R^,  is: 

R  :  Select  IT  iff 
3  i 


T  >  max  T  -  D,  j  <  -r  <  j  +  1 
2,1  '  l<r<k  J,r  '  L 


(5.1) 


and  D  is  chosen  to  satisfy 


inf  P{CS | R  >  -  P* 
*1 


(5.2) 


where  is  set  of  all  k-tuples  ^*^2’  satisfying  assumptions  given 

above . 

Now,  we  state  and  prove  a  theorem  related  to  the  infimum  of  probability  of 
a  correct  selection  when  rule  is  used.  Let  FU)  (x)  denote  the  distributions 

with  median  A^,  i  *  1,2,  ...,  k  . 

Theorem  5.1: 


If  F^j(x)  >  F^(x)  ,  for  all  x,  C(0)  =  4  and  G  ^^(x  +  A^) 
Is  nondoe re. is lng  (non  1  m: rens ing)  in  x  »  0  (x  <  0)  and 
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_d_ 

dx 


(x  +  A 


Ik] 


.  >  —  G(x)  |  >  0 

x^O  -  dx  'x=0 


then 


inf  r{CS|R  } 


/  gJ"-1  (t  +  D)  dGj(t) 


where  C,  is  defined  as  before. 
J 


Proof : 


By  stochastic  ordering  of  the  order  statistics,  we  have 


P{CS|R3>  >  /  Hk  *  (t  +  D)  dHj  k(x) 


(5.3) 


P|X.  ,  >  max  X,  -  D 


PNk  :  “i-'  ‘  Vl)'D| 


where  X,  ,  ,X.  ....  X.  .  are  i.i.d.r.v.  with  distribution  H,  ,  (x)  =  B  F.  , 

j,l  J « 2  j,k  j,k  j,n  [k] 

Let  *(x)  -  G_1F[k]  (x  +  A[k])  =  G”1  M(x)  where  M  =  BJ>n  F[k]  (x  +  &(k)) 

is  the  distribution  of  X.  -  A..  ,  .  Note  that  $(X.  -  Ari  .)  =  Y  has 

j,r  [k]  j  »r  [k]  j.r 

distribution  G,  .  Now  +  in  x  >  0,  +  in  x  <  0  and  1>,(0)  >  1 

j  x  -  x 

imply 


(l<r<k-l  (Xj<r  "  A[kj))  "  *(Xj»k  ~  A[k]) 
max  (X.  “  A . .  . )  —  (X,  i.  ~  A..  .) 


>  1 


(5.4) 


l<r<k-l 


j.r  u(k]'  Vrtj,k  u[k]' 


To  see  (5.4)  we  need  only  verify  that 


»(y)  -  bW  >  ! 
y  -  x 


for  y  >  x 


(x) 
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^  r  1. 1  (Af ,  i  ) 


Note  that  4>  *  (0)  -  L  _  *  where  f.^j  <g)  is  the  density  of  (G) 


Case  1: 


(0  <  x  <  y)  +  in  x  >  0  and  $'(0)  >  1 


imply 


$1*1  >  1 

y  -  x  -  x 


Case  2: 


(x  <  y  <  0)  follows  from  +  in  x  <  0  and  ^'(0)  >  1  . 


Case  3: 


(x  <  0  <  y)  .  In  this  case  $(x)  <  x  and  $(y)  >  y 


imply 


y  -  x 


We  have,  by  (5.4) 


max  Y.  -  Y.  .  >  max  X.  -  X.  . 
l<r<k-l  ^,r  J»r 


l<r<k-l 


(5.5) 


which  implies 


P  Y4  .  >  max  Y. 

I  3>k  -  X<r<k-1  3>r 


D  <  P  X.  .  >  max  X.  -  D 
-  3lk  l<r<k-l  3>r 


(5.6) 


which  proves  the  result. 


Note  that  we  need  not  have  restricted  attention  only  to  medians. 
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6.  SKl.F.CTION  WITH  RKSPFCT  TO  THE  MEANS 

Let  p^  be  the  mean  of  the  distribution  F(x;  p^),  i  =  1,2,  ....  k  and 
assume 


(a)  F (x;  Vjjj)  >  F(x;  y[k]^  for  1  =  1,2»  •••»  k  “  1 

(b)  F(x;  p^j)  <  G  (x)  =  1  -  e  X  for  i  =  1,2,  ...,  k  . 

Note  that  by  assumption  (b)  we  arc  confining  attention  to  the 
of  distributions.  It  will  also  be  convenient  to  assume  F(0; 
all  i  . 


and  all  x 


so-cal 1 ed 


IFR  class 
for 


n 

Let  x  =  £  x  /n  ,  where  x  is  the  j-th  observation  in  a  random  sample 

j=l  J  J 

of  size  n  from  .  Let  K^(x)  =  K(x;  p^)  be  the  distribution  of  x^  .  Then 
if  (x)  *=  K (x;  p^j) 

K[^j(x)  >  (x)  for  i  =  1,2,  ....  k  -  1  and  all  x  .  (6.1) 

K[i]  -<  G  for  i  -  1,2 . k  .  (6.2) 


(6.1)  is  an  immediate  consequence  of  (a)  while  (5.2)  follows  from  (b)  and  the 
closure  of  IFR  distributions  under  convolution  (see  Barlow,  Marshall  and  Proschan 


(1963)). 

If  we  are 
the  largest  p 


interested  in  selecting  a  subset  containing  the  population  with 
jkj  ,  we  use  the  rule 


:  Select  population  iff 

x  >  c '  max  x 

1  "  Ilk  1 
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It  foil ows  l ha t 

1  hoot  or  (< .  1  : 

or* 

i*U:s|k4»  :  /  [c;(-x-->  Jk_  1  dc:(x) 
o  c ' 

who  to  G(x)al-eX 

The  proof  is  the  same  as  for  Theorem  2.1.  The  disadvantage  is  that  the 
right-hand  side  of  the  inequality  is  independent  of  n  .  However,  by  restricting 
the  class  of  distributions  to  the  gamma  family  we  can  obtain  a  lower  bound  which 

depends  on  n  . 


App 1 i oat  i on  to  the  Selection  of  gamma  Populnt i ons 


Lot  us  consider  k  populations  with  densities 


1  x 

X®  xa_i  e  /  f (a) ,  x  >  0,  >  0,  i  -  1,2 . k 


Assume  that  a  >  1  ,  but  otherwise  unknown.  This  implies  that  the  distributions 
are  IKK.  Wo  are  interested  in  selecting  the  population  with  the  largest  (smallest) 
value,  '[k]  O(ij)  .  based  on  an  independent  sample  of  size  n  from  each  of 
the  k  populations.  Note  that  u  ^  j  =  a/A^j  for  i  =  1,2,  ...,  k  .  The  subset 

selection  rule  based  on  the  sample  means,  x . ,  i  =  1,2,  ...,  k  ,  is  R.  as  before. 

i  4 

(a)  r 

Let  Cv  denote  a  gamma  distribution  with  parameter  a  .  Since  l  x.. 

j  =  l  1J 

(nci) 

is  distributed  as  a  gamma  random  variable  with  distribution,  G v  it  follows 
from  a  result  of  van  Zwet  (1964)  that 


G(nCl)  -<  G(n) 


c 


vlii'ii  u  •  1  .  It  follows  lli.it  ill  this  case 


P{Cs!R4)  >  /  [G(n)  (— )Jk_l  dG(,0(x)  . 

o  o' 

The  constant  c'  is  determined  by 


(6.3) 


/  |G(n)(~)]k”1  dG(n)  (x)  ■=  P*  .  (6.4) 

o  c ' 

The  values  of  c*  are  tabulated  in  Gupta  (196'i).  It  should  he  point'd  out 
that  for  select  inp,  llie  population  with  the  smallest  A  ,  the  rule  can  le 
modified  to: 

:  Select  population  iff 

x  <  d  min  x, 

1  "  l<i<k  1 

where  d  is  determined  by 

J  [1  -  G<n)(*)]k_1  dG(n)(x)  *=  P*  . 
o 

The  values  ol  d  are  tabulated  in  Gupta  and  Soho)  (19b?). 

The  shape  parameter  u  *  l  need  not  be  the  same  for  all  populations.  It 
is  only  necessary  that  the  distribution  of  the  population,  n^j  •  wtth  the 
largest  mean  be  stochastically  larger  than  the  others. 
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